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An extended time-dependent Hartree-Fock theory, known as the time-dependent density-matrix 
theory (TDDM), is solved as a time-independent eigenvalue problem for low-lying 2 + states in 24 
to understand the foundation of the rather successful time-dependent approach. It is found that 
the calculated strength distribution of the 2 + states has physically reasonable behavior and that 
the strength function is practically positive definite though the non-hermitian hamiltonian matrix 
obtained from TDDM does not guarantee it. A relation to an extended RPA theory with hermiticity 
is also investigated. It is found that the density-matrix formalism is a good approximation to the 
hermitian extended RPA theory. 
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I. INTRODUCTION 



The time-dependent density-matrix theory (TDDM) Q is a version of extended time-dependent Hartree-Fock 
theories which dynamically include two-body correlations. Although TDDM has been formulated to deal with large 
amplitude collective motions 0, 01, it has also been applied to small amplitude oscillations: Dipole and quadrupole 
giant resonances in stable nuclei 0, Q and low-lying states in unstable oxygen isotopes 0, ■ Although the obtained 
results are quite reasonable, the foundation of TDDM for small amplitude motions is not clear in its time-dependent 
form. In this paper an eigenvalue problem of the Small amplitude limit of TDDM (STDDM) is solved for low-lying 2 + 
states in O to better understand the time-dependent approach. The hamiltonian matrix of STDDM is not hermitian 
as will be shown below. STDDM is compared with an extended RPA (ERPA) theory which has hermiticity || and its 
relation to the ERPA theory is clarified. This paper is organized as follows: The formulation of STDDM is presented 
in Sect. 2. A numerical solution of STDDM for low- lying 2 + states in 24 O is shown in Sect. 3. A comparison of 
STDDM with a more elaborate ERPA with hermiticity is made in Sect. 4 and Sect. 5 is devoted to a summary. 



II. SMALL AMPLITUDE LIMIT OF TDDM 



TDDM gives the time-evolution of the one-body density-matrix p(l, 1') and the correlated part C(12, 1'2') of the 
two-body density-matrix, where numbers denote space, spin, and isospin coordinates. The equations of motion for p 
and C have been obtained by truncating the well-known Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy of reduced 
density matrices 0- STDDM has been formulated according to the following steps 0: 1) Linearizing the equations 
of motion for p and C with respect bp and SC, where bp and SC denote deviations from the ground-state values po 
and Co i.e. Sp = p — po and SC — C — Co, respectively. 2) Expanding Sp and SC with single-particle states ip a as 

Sp(U',t) = ^w(t)i(l,i)e(l'*). (1) 

eta' 

<5C(121'2',i) - J2 ^« 1 «»«x(*)V'a 1 (l,t)V'a 2 (2 ) *)^ i (l , ,*)V'^(2 / ) *) ) (2) 
where ip a is chosen to be an eigenstate of the mean field hamiltonian /io(po) : 

h (p )iP a {l) = -^—i) a (l) + J d2v(l, 2)M2, 2)Vai)-p (l,2)Va(2)] =e Q Va(l)- ( 3 ) 
Here the one-body density matrix po is given as 

Po(H') = £nL^(l)C'(l')- (4) 

a a 1 
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After Fourier transforming the linearized equations, the equations for x aa i{uS) and X ai0l20l i a i (u>) can be written in 
matrix form, 



a c 
b d 



(5) 



where u> is an eigenvalue. The matrices a, b, c, and d are explicitly given in the Appendix. Eq.JSJ) can also be obtained 
from the following equations: 



{$o\[a+>a a ,H]\<i>) = uj($ \a+,a a \$>), 
($ \[a+, a+, a a2 a ai , H]\$) = uj($ \a+, a+, a a2 a ai \$), 



(6) 
(7) 



where [ ] stands for the commutation relation, H is the total hamiltonian consisting of the kinetic energy term 
and a two-body interaction, |$o) is the ground-state wavefunction and \$) the wavefunction for the excited state 
with excitation energy ui. Linearizing Eqs.© and J7J) with respect to x aa > — ($o| a~^,a a \$>) and A" Q ^ Q /^< = 
($o|a^' a a' a a 2 a a 1 \&) , and using the occupation matrix rfi aa , = {<&Q\cL^,a a \<&Q) and the two-body correlation ma- 
trix C° iQ2Q / Q ,/ = (^o\ a a' a a' aa 2 a ai\^o) — • / ^( Tl a 1 a' n a a'): where A is an antisymmetrization operator, we arrive at 
Eq.JSJ. Eq.JSJ is reduced to the second RPA (SRPA) |llLll2| under the two assumptions: The one is the Hartree-Fock 
(HF) approximation for the ground state, in which rP aa , — S aa >f a with f a = l (0) for occupied (unoccupied) states 
at zero temperature and C„ ia2a / Q < =0- The other is to restrict X aia2Ct ^ a ^ only to two particle - two hole components 
and their complex conjugates. 

The hamiltonian matrix of Eq.JSJ) is not hermitian, i.e. b ^ c + , as is easily understood from its explicit form. For a 
non-hermitian hamiltonian matrix, the left-hand-side eigenvectors of the hamiltonian matrix constitute a basis which 
is orthogonal to (x aa i, A Q , lQ2Q ' iQ ,^), and the ortho-normal condition is written as 



(8) 



where \fi) represents an eigenvector (x aa ', X axa2a ^ a i ) with the eigenvalue and |/2) the left-hand-side eigenvector 
of the hamiltonian matrix with the same eigenvalue: 



(£* X* 



a c 
b d 



uv(F X*). 



The completeness relation is written as 



where / is the unit matrix. In general the strength function is defined as 

S(E) = \^^\Q\^o)\ 2 S(E-E^, 

-B„>0 



(9) 



(10) 



(11) 



where \^q) is the ground-state, is an excited state with an excitation energy E^, and Q an excitation operator. 
The strength function in STDDM is calculated in a way similar to that used in TDDM. In TDDM, an initial excited 
wavefunction |$(t = 0)) is given by boosting the ground-state wavefunction |$o) as 

|$(t = 0)> =e lfe(3 |$o). (12) 
Then the initial values of the time-dependent amplitudes x aa i and X aia2a ' i0l ^ are given at first order of k as: 

W(0) - (^(t = 0)\a+a a mt = 0)}^ik{^o\K,a a ,Q}\^ ) (13) 
X aia2 a' ia ' 2 (0) w ik{$ \[a+a+a a2 a ail Q]\<5>Q). (14) 

In TDDM the time-dependence of x aa > and X ai „ 2a / a / is determined by solving numerically the TDDM equations. 
It can also be expressed using the eigenstates of Eq. JSJ as 



x(t) 
X(t) 



exp 



a c 
b d 



t/h 



i(0) 
X(0) 



= ^ exp 



.Ufj.t 



X„ 



x(0) 
X(Q) 



(15) 
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where the completeness relation Ea. (llO|) is inserted. The time-dependent moment Q(t) = ($(t)\Q\Q(t)) for a one-body 
operator becomes 

Q(t)= \ ex P ) lj2^ a ^ a '^ Xa ' a ^n (^ X (°) + K X (°)) + complex conjugate I , (16) 

fle(u^>0) I ^ ' \aa' J ) 

where we use the fact that both states with oj^ and — w* are eigenstates of Eq.©. The strength function in TDDM 
is obtained from the Fourier transformation of Q(t) as 



S{E) = ^- f° Q(i) e -rt/2R sin ^* dt 
irkn J Tl 



(17) 



where an artificial damping factor T is introduced to obtain a smooth distribution for S(E). Equations l|16fl and 117|) 
suggest that the strength function in STDDM has the form 



Re(uj u )>0 \aa' ) \ BB' 



1 



E-UJ^+iT/2 



J2(u\Q\a')x a 'a(v) ) I YsWWH'P^) 1 E + u*+iT/2 ]} > (18) 



where x l aa , (p) is defined as 



5L'(M) =5Z^°l[ a a' a «' a A a A']|*o)iAA'(/i)+ J! (^olIoiacO^O^OASOA'JI^O^AaAi^O*). (19) 
AA' AiA2A^A^ 

The strength function S(E) in STDDM is not guaranteed to be positive definite, as is easily understood from its 
expression Eci. (|18fl . 

III. NUMERICAL SOLUTION 

In this section we present a numerical solution of Eq.JSJ calculated for low-lying 2 + states in 24 0, where a time- 
dependent approach has previously been applied • The E2 strength function is calculated according to the following 
steps: 

1) A static HF calculation is performed to obtain the initial ground state. The neutron 2s\/2 state is assumed to 
be the last fully occupied neutron orbit of 24 O. The Skyrme III is used as the effective interaction. It has been used 
as one of standard parameterizations of the Skyrme force in nuclear structure calculations even for very neutron rich 
nuclei [lit Il4j. The single-particle wavefunctions are confined to a sphere with radius 12 fm. The mesh size used is 
0.1 fm. 

2) To obtain the correlated ground state l^o): we evolve the HF ground state using the TDDM equations for 
ip a (l,t), n aa i(t) and C aia2a ' i0l ' 2 (t) 0,Q| and the following time-dependent residual interaction 

v(t) = (1 - e-'/Xr - r'). (20) 

The time constant r should be sufficiently large to obtain a nearly stationary solution of the TDDM equations • 
We choose r to be 300 fm/c. In the calculation of n aa i(t) and C ai0l2a ' ia > 2 (t), the single-particle states are limited to 
the neutron orbits 1(^/2, 2s\/2 and 1^3/2- In a consistent calculation the residual interaction should be the same as 
that used to generate the mean field. However, a Skyrme-type force contains momentum dependent terms, which 
make the computation time of two-body matrix elements quite large. Therefore, we need to use a simple force of the 
5 function form v oc 6 3 (r — r'). We use the following pairing- type residual interaction of the density-dependent 5 
function form 0] 

v(r - r') = v (l - p(r)/p Q )S 3 (r - r'), (21) 

where p(r) is the nuclear density. The parameters po and vo are set to be 0.16fm~ 3 and —1000 MeV fm 3 , respectively. 
Similar values of po and i>o have been used in the Hartree-Fock-Bogoliubov calculations 0, 0, ^| in a truncated 
single-particle space. The time step used to solve the TDDM equations is 0.75 fm/c. 
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FIG. 1: Strength distributions of the neutron quadrupole modes in 
and SRPA (thin dotte line). 



'O calculated in STDDM (solid line), RPA (dotted line) 



3) At t = 5r we stop the TDDM calculation and solve Eq.JSJ. Since we are interested in low- lying states which 
originate in inner-shell transitions, we use only the the neutron lds/2, 2si/2 and 1^/2 states. The dimension of the 
hamiltonian matrix is about 650 x 650 when these single-particle states are used. The strength function for 2 + states 
is calculated using Eq.^HJ with Q = r 2 Y- 



20 1 



We use T = 0.5MeV. 



The eigenvalues of some 2 + states become imaginary because the hamiltonian matrix of Eq.© is not hermi- 
tian. However, their imaginary components are quite small (less than 0.08MeV). Some 2 + states have also negative 
quadrupole strengths because positivity of S(E) is not guaranteed. However, the negative contributions are so small 
that S(E) becomes almost positive in the entire energy region when it is smoothed with T — 0.5MeV. The obtained 
result in STDDM (solid line) is shown in Fig.l, where the strength functions in RPA (thick dotted line) and SRPA 
(thin dotted line) are also presented for comparison. The first 2 + state calculated in STDDM is energetically shifted 
upward and becomes significantly more collective as compared with that in RPA. The increase in the excitation energy 
is due to the lowering of the ground state, and the enhancement of the collectivity of the first 2 + state is due to the 
mixing of two-body configurations induced by ground-state correlations. These properties of the first 2 + state under 
the influence of ground-state correlations are similar to those obtained from quasi-particle RPA calculations [20I l2l) 
and also from the TDDM calculation 0- Large single-particle space including continuum states has been used in 
these realistic calculations H3, EI so as to deal with giant resonances as well as low-lying states. Therefore, only 
a semi-quantitative comparison between our previous TDDM result and the present STDDM calculation, which 
can only be performed in very truncated single-particle space, should be made for lowest-lying states. In SRPA where 
ground-state correlations are neglected, the 2 + states become also collective due to the coupling to two particle - two 
hole configurations, as shown in Fig.l. However, their excitation energies are simply shifted downward. This might 
cause serious problems when unperturbed particle - hole energies are small and (or) when the coupling to two-body 
configurations is strong. The numerical solution of Eq.(|5jl shown in this section is physically reasonable and also gives 
additional justification to its equivalent time-dependent approach. 



IV. DISCUSSIONS 



In this section we discuss a relation between STDDM and a more elaborate extended RPA theory (ERPA) which 
has hermiticity |§J. First we briefly present the ERPA theory. The ground state |$o) m ERPA is constructed so that 



($ Q \[H,a+a a ,]\$ o )=0, 



(22) 



($ |[i?,a+ 0+ 0^0^] |$ ) = 0, 



(23) 



($o|[-H>+ a+ 2 a+ s a a ja a £<vJ|$o) = 



(24) 
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are satisfied for any single-particle indices. In other words n° ,, C° / , and the three-body correlation matrix 

J o L OLOL ' aia2Q 1 Q 2 J 

^aia 2 ct3a' q' a' are determined so that the above equations are satisfied. The explicit expressions for Eas. l|22|l - I|24|) are 
shown in Refj^]. Although it is not evident to find an analytic solution of Eas. l|22[l - ll24|) 22], a method for obtaining 
nj , and C° 1Q2Q ' a < numerically has been proposed in Ref.[n| and already been tested for realistic nuclei [UOH^ as 
was shown in Sect. 3. The ERPA equations are formulated using the equation of motion approach [24l l25j as 

(H> \[[aia a ,,H},Q + }\y Q ) = uj(y Q \[a+a a ,,Q+}\y Q ) (25) 



<*o|[[: «a a >a a[ :,H], Q+]\9 ) = w(9 \[: a+a+a^, :,Q + ]|*o), (26) 
where the operator Q + is defined as 

Q + = J2( Xxx ' a * ax ' + X AiA 2 a;a^ : a+^a^ax^ :), (27) 

and l^o) is assumed to have the following properties 

Q+|vf ) = |*) (28) 
Q|*o> - 0. (29) 

In Eqs.@ and : : stands for : a^a^a^a^ := a^a^a^a^ - A(a+ 1 a a ! L ($o\ai 2 aa' a \®o) + 

a a 2 a <x' 2 (^o\ a ai a a' l^o))- Equations l|25|) and (|26[) can be written in matrix form 



(30) 



where each matrix element is evaluated using the ground state |$o) 

5i(a'a:AA') = (Q \[a+a a >, 4 a A']|$o), (31) 

S 2 (a' 1 a 2 a 1 a 2 : A^A'^) = ($ |[: a+a+ 2 a Q ^a ai :, : a^a^a A ^a A; :]|$ >, (32) 

Tx(a'« : A1A2A1A2) = ($ \[a+a a ,, : 0+0+°^°*! :]|*o>, (33) 

r 2 (aia 2 aia 2 : AA') = ($ |[: a+ a+ 0^0^ :, a^av]|$ ), (34) 

A(a'a :AA') = <*o|[[a+a«/, fl],aJaA']|$o), (35) 

B(a' 1 a' 2 aia 2 : AA') = ($ |[[: a+a+a^flai :, #], a+a v ]|$ ), (36) 

C(a'a : A 1 A 2 A' 1 A 2 ) = <$ |[[4<W, ff], : a+ a+ a A ^a Vi :]|$ ), (37) 

£)(a' 1 a 2 aia 2 : AiA 2 A' 1 A' 2 ) = ($ |[[: a+a+aQ^Q^ :, i/], : a+ a+ a^a^ :]|$ }- (38) 

All matrices in the above are written in terms of n° , , C° / / and C° > > t , which are shown in Rcf. 81. Due 
to Eas. (|22ll - 124l) . the above matrices have the following symmetries 

A(a'a : AA') = A(X'X : aa') = A(W : a'a)*, (39) 
B(a' 1 a' 2 aia2 : AA') = C(A'A : aia 2 Q;ia 2 ) = C(AA' : aia 2 aiQ; 2 )*, (40) 
Z?(a' 1 a 2 aia 2 : AiA 2 A^A 2 ) = D(X' 1 X' 2 XiX 2 : a.\a 2 ol^OL^) = D(X\ X 2 X[ A 2 : a' 1 Q! 2 aiQ; 2 )*. (41) 

We explain in more detail the last relation as an example. The following relation always holds due to the operator 
identity 

D{a' 1 a' 2 aia 2 : AiA 2 A' 1 A 2 ) - D(X' 1 X' 2 X 1 X 2 : aia 2 aia 2 ) = ($o|[[: a+ a+ a^ae,/ :,#],: a^a^a^a^ :]|$ ) 

- <*o|[[: a^a^o^aA'j :,#], : a+ a+o^Oc,/ :]|$ ) 
= (*o|[#, [: «fla;a„; :, : 0^,0^0^ :]]|$o). (42) 

The commutation relation between two-body operators in the last line of the above equation is reduced to three-body 
operators at most. Since Eas. l22|) - (|24|l are valid, the last line of Ea. (|42|) is identical to zero. Thus, the symmetry of 
the matrix D is proven. 
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The ortho-normality condition in ERPA is given by [27 

Si T\ 



K>K>)\ r s 2 (x )=V. ( 4 3) 



where a; M and A M constitute eigenstates of Eq. 1)30(1 with w =w M . The completeness relation becomes 



The transition amplitudes for one-body and two-body operators, z aa > = ('$>o\a~^,a a \' ( &) and Z aia ^ a i^ = (^>o 
a a' a a' a ai a a 2 '■ l^ 7 ): are calculated as follows 



z \ _ ( Si Ti \ ( x 
Z -\T 2 S 2 { X 



(45) 



Now we show a relation between STDDM and ERPA. When the eigenvector (x, X) in STDDM is transformed to 
(y, Y) as 



* ^ _ f Si Ti\(y 

x -\t 2 S' 2 Y 



(46) 



where S 2 has no terms with the three-body correlation matrices, the equation in STDDM becomes 

(47) 



aSi +cT 2 aTi + cS^\ ( y \ _ ( S x T x \ ( y 
bSi + dT 2 bTi + dS 2 ) \ Y ) u \T 2 S' 2 ) \ Y 

It is straightforward to prove that when the three-body correlation matrices C aiot2a3Ct i a ' 2a i are neglected, three blocks 
of the above hamiltonian matrix satisfy the identities A = aSi + cT 2 , B = bSi + dT 2 , and C = aT\ + cS 2 . We 
explain the relation B = bSi + dT 2 in detail. When two-body and three-body operators are decomposed using 
a tt 1 a tt 2 a a^ a a' 1 = : a+ a+ a a ^a Q ^ : -M(a+a a / ($ |a+ |$ ) + a+ (^ok+a^ |$ )) for a two-body operator and a 
similar relation for a three-body operator, the commutation relation [: a +, a a > a a ' :,H] becomes 



H] = 2_,^{ a i a 20i'ia' 2 : W')a\,a\ 



*ai"a2 a 2 a i '' J / j 

AA' 

+ 2_j d(aia 2 a[a 2 : A1A2A1A2) : ay^ay^a\ 2 a\ ± : 

Ai A 2 A.^ A 2 

+ £ e{aia 2 a' 1 a' 2 : Ai A2A3A1A2A3) : ay^a y ^aya\ s a\ 2 a\ 1 :, (48) 

AiA2A3A 1 A 2 Ag 

where the matrix e consists of matrix elements of the two-body interaction. In STDDM, the last term on the right- 
hand side of the above equation is neglected. It is clear from Eq. that B = ($o|[[ : a a 1 a a 2 a a' 2 a a' 1 : , H], a x a\<]\$o) = 
bSi+dT 2 . The matrix D has some terms which cannot be written by using bT\+dS 2 even if the three-body correlation 
matrices are neglected: For example, all terms with the square of the two-body correlation matrix cannot be expressed 
with bT\. Thus it is found that STDDM is identical to ERPA under the following two assumptions: The first is that 
the three-body correlation matrices can be neglected and the second is that D « bT\+dS 2 . Our numerical calculations 
so far performed suggest that these assumptions are quite reasonable. 



V. SUMMARY 



The eigenstates of STDDM were calculated for the low-lying 2 + states in O to better understand the foundation 
of its corresponding time-dependent approach. It is found that STDDM properly deals with the effects of ground-state 
correlations on the low-lying 2 + states and that the non-hermiticity of STDDM is quite moderate: The eigenvalues 
have quite small imaginary parts and the strength function is practically positive definite although it is not guaranteed 
because of its non-hermitian form. A comparison of STDDM and an extended RPA (ERPA) theory with hermiticity 
was also made, and it is found that STDDM is a reasonable approximation to the ERPA theory. 



APPENDIX A 

The matrices in Eq.JSJ) are given by: 

a(aa' : AA') = (e Q - e a ,)S aX S a , x , - ^((/3\'\v\a'\) A n° a p - (a\'\v\p\) A n%,), (Al) 



b(a 1 a 2 a' 1 a' 2 : AA') = S ai x(^2[(5 a20 - n^)n° a ,^ a , + n° 2(3 (<5 7Q < -n° a; )(<W 2 - ™L 2 )] (A'/^^U 

+ E(( A ' a 2bl/37>^ 7Q , Q , + (A'/3|«|ai7)A^ 27Q ^ - (X' P\v\a' l7 ) A C o a2ja , } 

Pi 

+ W£[(*ai/»-<,j)n^ng a , +n|^ 

Pit 

+ Y^ti X ' a Mfr}C o 07a , a , + (X'PIvK^aC*^ - (X'P\v\a' 2l ) A C° ai7a ,J} 

Pi 

+ 8 a > iX >{Y^K S Sa> 2 ~ n° Sa , 2 )n a ail3 nl 2j + n£ Q , {S aiP - n° ai/3 )(S a2l - n° 27 )] ((3j\v\\6) A 

PlS 

+ ^[(/?7b|A4)C° lQ2/37 + (ai/3|«|A 7 )AC*« 27a , /3 - {a 2 p\v\X 7 ) A C° aiW J} 

Pi 

~ 8 a2 \>Q2l( S S<x'i ~ n L' 1 ) n " 1 /3 n a 27 + ~ n ° ai p)( S c<2l ~ n° a2l )]{/3-/\v\XS) A 

Pl$ 

+ ]T[</?7MAai}C° lQ2/97 + (aM^)AC° a2Wl p ~ ^v^) A Cl il<x ,J} 
Pi 

+ EK«i^>I/3A)aC° Q2Q , q , - (a 2 \'\v\p\) A C° 0aia , a , 

P 

(P\'\v\a' 2 \) A C° aia2a , + (PX'Ma'iXUCl^J, (A2) 
c(aa' : A 1 A 2 A' 1 A^) = (a^|«|A 1 A a )* a 'Ai - <A' 1 A>|a'Aa}*aA 1) (A3) 
d(aia 2 Q;' 1 Q;2 : \ 1 \2\' 1 \' 2 ) = (e ai + £a 2 _ - ea 2 )^aiAi^a 2 A 2 ^aiAi^a 2 A 2 

/3< 27 ){/37l«|AiA 2 ) 

/3 7 
187 

+ S Q2 a 2 5 Q2 a 2 53((«iAi|w|/3Ai) A ^ Q , - OSAilvlaiAi^n"^) 
P 

+ S a2 x 2 S a[X[ ^((aiA 2 |i;|/3Ai) A ^ Q , - (/3A 2 |w|a 2 Ai) A < 1/3 ) 

/3 

+ ^aiAi^Ai £((a2A 2 M/3A 2 } A r^ Q , ~ (/ 3A 2l w l"2 A 2)A< 2/3 ) 

+ <WA 2 A 2 £«a 2 AiM/?A 2 ) A n° Q < - (/3A' 1 |t;|«' 1 A 2 ) A < 2/3 ), (A4) 

/3 

where n° Q , = (^olaiaal^o) and C aiOi2a ,^ = ($o|a^a+, a^a^ |$ ) - •^( n a 1 a' 1 n a 2 a !2 )' and the subscript ^ indicates 
that the corresponding matrix element is antisymmetrized. 



[1] M. Gong and M. Tohyama, Z. Phys. A335, 153 (1990). 

[2] M. Gong, M. Tohyama and J. Randrup, Z. Phys. A335, 331 (1990). 

[3] M. Tohyama and A. S. Umar, Phys. Rev. C65, 037601 (2002). 

[4] F. De Blasio et al., Phys. Rev. Lett. 68, 1663 (1992). 



8 



[5] M. Tohyama, Nucl. Phys. A657, 343 (1999). 

[6] M. Tohyama and A. S. Umar, Phys. Lett. B516, 415 (2001). 

[7] M. Tohyama and A. S. Umar, Phys. Lett. B549, 72 (2002). 

[8] M. Tohyama and P. Schuck, in preparation. 

[9] S. J. Wang and W. Cassing, Ann. Phys. 159, 328 (1985); W. Cassing and S. J. Wang, Z. Phys. A328, 423 (1987). 

[10] M. Tohyama and M. Gong, Z. Phys. A332, 269 (1989). 

[11] J. Sawicki, Phys. Rev. 126, 2231 (1962); J. Da Providencia, Nucl. Phys. 61, 87 (1965). 

[12] C. Yannouleas, Phys. Rev. C35, 1159 (1987). S. Drozdz et al., Phys. Rep. 197, 1 (1990). 

[13] T. Otsuka, N. Fukunishi and H. Sagawa, Phys. Rev. Lett. 70, 1385 (1993). 

[14] E. Khan and N. Van Giai, Phys. Lett. B472, 253 (2000). 

[15] M. Tohyama, Prog. Theor. Phys. 94, 147 (1995). 

[16] R. R. Chasman, Phys. Rev. C14, 1935 (1976). 

[17] J. Terasaki, H. Flocard, P. -H. Heenen, P. Bonche, Nucl. Phys. A621, 706 (1997). 

[18] T. Duguet, P. Bonche and P.-H. Heenen, Nucl. Phys. A679, 427 (2001). 

[19] M. Yamagami, K. Matsuyanagi and M. Matsuo, Nucl. Phys. A693, 579 (2001). 

[20] M. Matsuo, Nucl. Phys. A696, 371 (2001). 

[21] E. Khan, N. Sandulescu, M. Grasso and N. Van Giai, Phys. Rev. C66, 024309 (2002). 

[22] M. Tohyama, P. Schuck and S. J. Wang, Z. Phys. A339, 341 (1991). 

[23] M. Tohyama, Phys. Rev. C58, 2603 (1998). 

[24] D. J. Rowe, Rev. Mod. Phys. 40, 153 (1968). 

[25] P. Ring and P. Schuck, The nuclear many-body problem, (Springer, 1980). 

[26] M. Tohyama, Z. Phys. A335, 413 (1990). 

[27] K. Takayanagi, K. Shimizu and A. Arima, Nucl. Phys. A477, 205 (1988). 



